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Abstract 

We calculate the virtual and bremsstrahlung fermionic corrections due to the interfer- 
ence of the electro- and chromomagnetic dipole operators in the inclusive B — > Xgj 
decay at O(a^) ^^'^ present analytical results for both the total decay rate and the 
photon energy spectrum. 



1 Introduction 



The present experimental world average of the branching ratio of i? — > X^^y, which includes 
measurements by CLEO, BaBar and Belle [1-3], is performed by the Heavy Flavor Averaging 
Group [4] and, for photon energies > 1.6 GeV, is given by 

Br{B X,7) = (3.52 ± 0.23 ± 0.09) x 10"^ , (1.1) 

where the errors are combined statistical and systematic and due to the extrapolation to the 
common lower-cut in the photon energy, respectively. Moreover, the total uncertainty, being 
already below 7%, is expected to reduce down to 5% at the end of the B-factory era. 

In order to keep pace with the improving experimental accuracy the theoretical prediction 
of the B Xs7 branching ratio has to be known at the next-to-next-to-leading order 
(NNLO) level. A first estimate of the B — ^ Xg'y branching ratio at this level of accuracy has 
been presented in [5]. For Ey > 1.6 GeV it reads 

By{B X,7) = (3.15 ± 0.23) x 10"^ . (1.2) 

This estimate includes the three- loop dipole operator matching conditions [6] , the three-loop 
mixing of the four-quark operators [7], the three- loop mixing of the dipole operators [8], and 
the four-loop mixing of the four-quark operators into the dipole operators [9] . Also the two- 
loop matrix elements of the electromagnetic dipole operator together with the corresponding 
bremsstrahlung terms (at rric = 0) [10-13], as well as the three-loop matrix elements of the 
four-quark operators within the so-called large-/?o approximation [14] have been taken into 
account. Finally, in order to obtain estimates of these matrix elements (together with other 
still unknown ones) at the physical value of the charm quark mass rric an interpolation in 
the latter has been performed in [15]. 

We should mention here that there are several perturbative and non-perturbative effects 
that have not been considered when deriving the estimate given in (11.21) . Some of them 
are already available in the literature: the four-loop mixing of Oi, . . . ,0q into Os [9]; the 
bremsstrahlung contributions of the (02,02), (02,07) and (O7, 08)-interference at 0(a^/3o) 
[16]; the exact charm quark mass dependence of the (O7, 07)-interference at O(a^) [17]; the 
three-loop virtual corrections due to charm and bottom quark loop insertions into gluon 
propagators in the (Oi,07) and (O2, 07)-interference [18]; the updated knowledge of the 
semileptonic normalization factor [19-21]; photon energy cut-off related effects [22-25]; and 
estimates for the 0(asAQCD/'"^b) corrections [26]. Other effects are unknown at the moment, 
like the complete virtual and bremsstrahlung contributions to the (O7, Og) and (Og, Og)- 
interference at O(a^) (only the contribution of the (O7, Og) -interference at 0(a^/5o) is known 
[14, 16]), and of course the exact mc-dependence of various matrix elements beyond the 
large-/5o approximation, in order to improve (or even remove) the uncertainty due to the 
interpolation in [15]. The individual contributions listed above are all expected to remain 
within the uncertainty given in (11.21) . nevertheless they should be taken into account in future 
updates. 

In the present paper we repeat the calculation of the (O7, Og)-interference contribution 
performed in [14, 16] and extend it to include not only the effects of massless quark loops 
but also those due to massive ones. More precisely, we calculate those 0{a1) contributions 
which can be obtained from the Feynman diagrams contributing to the (O7, Og)-interference 
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at 0{as) when dressing the gluon propagators with massless up, down and strange quark 
loops as well as with massive charm and bottom quark loops. We work out the effects of 
these contributions to the photon energy spectrum dr{b —>■ Xf^^^°^^'^'-f)/dEy and to the total 
decay width r(6 —>■ X^'^'^^°'^^'^'j)\e^>Eoj where Eq denotes the cut in the photon energy. 

The organization of this paper is as follows. In section 2 we present our final results for 
the total decay width and the photon energy spectrum and describe briefly the calculation 
of the relevant Feynman diagrams. The numerical impact of the fermionic corrections on 
Br[B — > Xgj) is estimated in section 3. Finally, we summarize in section 4. 



2 Fermionic corrections 

Within the low-energy effective theory the partonic b X^j decay rate can be written as 

T{b - xr°^7)i..>i.o = ^^^-g^)"^' \v,,v;f Y.cf{^^)cf{^^)G.,{E,,^^) , (2.1) 

where rrib and rribifJ.) denote the pole and the running MS mass of the b quark, respectively, 
Cj^^(/i) the effective Wilson coefficients at the low-energy scale, and £"0 the energy cut in the 
photon spectrum^ 

As already anticipated in the introduction, we will focus on the function Gjs{Eq, fi) 
corresponding to the interference of the electro- and chromomagnetic dipole operators 

^7 = TTT^ (sa'^'PRb) E^, (2.2) 

and 

Os = {-sa''^PnT%) G% , (2.3) 

respectively. In NNLO approximation this function can be decomposed as follows, 

GM,f^) = ^Cpy«(^o,/i) + Cpr(^)(^o,/i) + 0(af) , (2.4) 
where 

F(2)(zo,/i) = GpY^'^^^\zo,fi) + GAY^'^''-\zo,fi) 

+ TnNLY^''^^\zo, fi) + rRiVHF(''""H^o, ^) + TrNvY^^^^''\zo, fi) . (2.5) 

Here, zq = 2£'o/m6, Nl, Nh and Ny denote the number of light (m^ = 0), heavy (m^ = rrib), 
and purely virtual (m^ = rric) quark flavors, respectively, as(/i) is the running coupling 
constant in the MS scheme, and Cp, Ga and Tji are the color factors with numerical values 
given by 4/3, 3 and 1/2, respectively. In this paper we present results for the functions 
with i = NL, NH, NV. The calculation of the functions Y'^^'^^Zq, fx) with i = CF, 
CA is the subject of another publication [27]. 

In this paper we assume that the products Cf{n) Cf{fi) are real. That is our formulas are not applicable 
to physics scenarios beyond the standard model which produce complex short distance couplings. 
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Figure 1: Two sample b quark selfenergy diagrams which are proportional to the number of 
light, heavy and purely virtual quark flavors and whose 2-, 3- and 4-particle cuts contribute 
to the b ^ S'-y (blue dashed lines), b — > S'-yg (orange dashed lines), and b S'jqq (green 
dashed line) transitions at 0{a1). See text for more details. 

The functions F^^'*) with i = NL, NH, NV appearing in ( 12. Sp receive contributions from 
the b — >• S7, b — > S'jg, and b S'jqq (g G {u,d,s}, = 0) transitions. The latter 
are contained in the b quark selfenergies which arise from those at 0{as) when dressing the 
gluon propagators with massless and massive quark loops|§ Two sample b quark self-energies 
containing cuts with two, three and four particles in the intermediate state are displayed in 
figure [H As far as the diagrams containing massive quarks in the fermion loop are concerned, 
like, e.g., the one given on the left-hand side of figure [H, we do not have to calculate cuts 
with four particles in the intermediate state since such cuts would run through the bottom 
or charm quark loop and (i) b quarks are of course kinematically not allowed to appear in 
the final state and (ii) events involving charmed hadrons in the final state are not included 
on the experimental side. On the other hand, for the diagrams containing massless quarks in 
the fermion loop like, e.g., the one given on the right-hand side of figure [H, the contributions 
from qq production (g G {u, d, s}), that is four-particle cuts running through massless quark 
loops, have to be taken into account. 

We work in d = A — 2e space-time dimensions to regularize ultraviolet, infrared and 
coUinear singularities, and adopt the renormalization prescription from [12,17]. Most of 
the renormalization constants necessary to render our results ultraviolet finite can be found 
there. The only exceptions are those which describe the selfmixing of Os at one-loop and the 
mixing of Os into O7 up to two- loops; they can be extracted from [28]. For some technical 
details concerning the evaluation of the two-loop integrals involving both the bottom and 
charm quark mass mi, and mc, respectively, we refer the reader to the end of this section. 

In order to obtain a compact presentation of our findings we split the functions 
with i = NL, NH, NV into two parts, namely 

y(2,0(^^,^) = y(2.)(o,^)-5y(2.)(^^,^), (2.6) 

where the first terms give always the contribution to the full inclusive decay rate, and the 
second ones correct for the fact that in the experiments a lower cut in the photon energy is 
applied. Performing the same splitting for the function Y^^'^ appearing in (12.41) . our findings 

^ There are also b quark selfenergies where the photon runs from Oj to the quark loop and which cannot 
be obtained from those at 0{as). However, it is easy to show that these diagrams cancel amongst themselves 
without performing loop and phase-space integrations (see Furry 's theorem). 



for the two individual contributions at 0{as) read 

yW(0,/.) = ^(29 - 27r^)+^L,, (2.7) 

5Y^'\zo, f^) = lzo {z'o + 24) - ^(^o - 1) ln(l - ^o) - ^ U^izo) , (2.8) 
while those at 0(0^) are given by 

Y^'''''^\0, ^^)^-f^ (328 - 13vr^) - | (l8 - tt^) - y + y Ca , (2-9) 

64 



^^^''^"HO,/^) = ^ (244 - 27V37r - 61n') - ^ (l8 - n') L, 



- y - I C3 + 32^3 01.(1), (2.10) 

r(^'-)(0, /.) = -i^ [157 - 279p - TT^ (5 + 9p' - 42^^/^)] - | (l8 - tt^) L, 

- ^ + 1^ (22 - vr^ + lOp) Inp + (8 + 9p^) In^ 

- i 1^'/^ - ^ (1 - - ^(19 - 46p) $2(p) 

32 64 32 

- -(13 + 14p) $3(p) - Y ^^(P^ - Y ^^^^^^2(1 - p) 

+ I (5 + 9/ + 14pV^) Li,(l - p) - il^//^ Li, (1 - Vp) 

+ ^Li3(l-p) + ^Li3(l-i), (2.11) 

(5y(^''^"H^o, fi) = Zo {7zl - 17zo + 238) - ^ 6Y^'\zo, p)L^ 

g 

+ ^ (-20 - Qzl + 80^0 - 75 + 67r^) ln(l - zo) 

16 16 

- —{zo - 1) ln2(l -Zo)-Y In'^o ln'(l - -^o) 

32 32 

- ^(32^0 - 8) Li2(2;o) - y ln(l - zq) U2{zo) 

+ y Li3(zo) - y Li3(l - zo) + y C3 , (2.12) 

5yM(zo,/^) = -^5y«(^o,/^)i^M, (2-13) 

SY^''-'^\zo, f^) = 5y«(zo,p)(2L^ - Inp) . (2.14) 
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In writing these equations we introduced the short-hand notations 



^2 



p=^ and L„ = lnf-^V (2.15) 



The definitions of the auxihary functions $n(p) as well as those of the polylogarithms Li„(2;) 
and the Clausen function Cl2{z) can be found in appendix [Bl The numerical value of the 
Clausen function at 2; = 7r/3 is approximately given by 1.014942, and Cs ~ 1.202057 is equal 
to Riemann's theta functions ({n) at n = 3. Equations (12.111) and (12.141) hold for p > 0. 

Turning now to our findings for the photon energy spectrum, we rewrite the function 
G-rsiEo, fj,) as an integral over the (rescaled) photon energy, 

G7s{Eo,p)= dz , z= '-. (2.16) 

In NNLO approximation the integrand can be written as follows, 

= ^c^ya)(,,^) + ( c^y(^)(.,p) + o{al) , (2.17) 

dz 4:71 \ 4:71 J 

where, in analogy to (12.51) . 

Y^'\z,fi) = TnNLY^''^'^\z,f,) + TrNhY^'^''''\z, p) + TnNyY^^^^''\z, p) + . . . , (2.18) 

with the ellipses denoting terms which are proportional to the colorfactors Cp and Ca- The 
next-to-leading order (NLO) function Y'^^\z^jj) is given by 

yW(z,p) = ril.(0,^)5(l - z) + ^6Y^'\z,p) , (2.19) 

dz 

where 6Y^^\z,fi) can be obtained from (12. 8p by replacing zq by z, and Y^J^^^{0, n) can be 
found in appendix |X1 The latter function summarizes the contribution of all 2-particle cuts 
entering the functions Y^^^zq, fx). The terms proportional to Nl, Nh and Ny appearing in 
(12.181) can be written in complete analogy to (I2.19p . 

y(2'^)(^,/i) = r£i(o,/i) 5(1 - z) + ^5y(2'^)(^,/i) , (2.20) 

with dY^'^'^^{z, n) given in (I2.12I) - (I2.14I) and ^2-cuts(05 appendix |Al Since the contribu- 
tions of the 2-particle cuts are by themselves free of infrared and collinear singularities it 
was not necessary to introduce plus-distributions in (I2.20p and (12.191) . 

We remark that the terms proportional to N^, that is the functions Y^'^''^^\0, fi) and 
6Y^'^'^^\zo, p), are already known in the literature [14,16] and we completely agree with the 
results given there. The functions Y^^^'\zo,n) with i = NH, NV are however newl 

In the remainder of this section we will summarize the technical details of the calculation. 
However, we refrain from repeating the algebraic reduction procedure of the 2-, 3- and 4- 
particle cuts of the three-loop b quark selfenergies to a set of so-called master integrals as 



would like to thank Christoph Greub for checking equation (jA.5p . which contains the contribution of 
the 2-particle cuts being proportional to N^, numerically. 
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well as from discussing appropriate parametrizations of the phase-space integrals here since 
this has already been done in great detail in [12,13]. Instead, we will briefly describe how we 
solved the non-trivial two-loop integrals involving the two mass scales rrib and m^. First, we 
introduced Feynman parameters in the standard way and performed the loop-integrations. 
Subsequently, we applied the Melhn-Barnes technique [29,30] based on the relation 

T{-s)T{\ + s), (2.21) 



(x + y)^ r(A) J 2m 

c 

where the integration contour C runs from —ioo to +ioo such that it separates the poles 
generated by the two F functions. In this way all powers of a sum of several terms could 
be replaced by one- or two-fold Mellin-Barnes integrals which made the integration over 
the Feynman parameters trivial. Finally, we closed the integration contours C sidewards 
by a half-circle with infinite radius and summed up the enclosed residues. In our case all 
infinite sums involving the mass ratio rric/ mb could be reduced to the inverse binomial sums 
given in [31], and we obtained solutions for all two-loop integrals which are valid for arbitrary 
values of rric/mb. We checked our analytical results for the master integrals for several values 
of nic/mb by numerically integrating over the Feynman parameter representations. 

Two other checks of our calculation are provided by taking the limits — > and 
rric —>■ mb of the contributions of the 2-particle cuts proportional to A'^v'; 



limrSrH0,/^) = rS^(0,/.) 

.(2,NV)/p, \ _ ■^(2,nh), 



n:z' (0, /i) = (o, z^) , (2.22) 

which reproduce our results proportional to A^^^ and Nh- Also the limit p — 1 of the 
complete expression Y^'^''^\zo, fi) for i = NV reduces to that for i = NH. We note, however, 
that it is not possible to take the limit p — of the complete expression for i = NV since 
we excluded the contributions with massive cc-pairs in the final state, and hence some Inp 
terms being present in Y^'^'^^'>{zo, fj.) parametrize infrared and collinear singularities. The 
last check concerns the asymptotic behavior for rric ^ In this limit our result for the 
complete expression with i = NV reduces to 



which is in agreement with the asymptotic form found in [6] (see equation (5.10) of that 
reference) . 



3 Numerical impact 

In the left frame of figure [2] we show the dependence of the functions Y^^^^^{0,mb) on the 
mass ratio p = ml/ml for i = NL, NH, NV. As can be seen, the function for i = NH (blue 
dashed line) differs from that for i = NL (red solid line) by a factor of about —0.5. On the 
other hand, at the physical value p ~ (0.262)^, the function for i = NV (black dotted curve) 
has a smaller value of about 15% compared to that for i = NL. 
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Figure 2: 1^2-cuL(05 '"^ft) f-'eftj and F'^^'*^(zo, '"^ft) (right) as a function of p = ml/ml for i = NL 
(red solid line), i = NH (blue dashed line) and i = NV (black dotted curve). The vertical 
lines indicates the physical value p ~ (0.262)^. 



The right frame of figure [2] displays the functions Y^'^''^\zQ,mb) for i = NL, NH, NV 
as functions of the mass ratio p. They differ from the ones shown in the left frame by 
adding the contributions from the 3- and 4-particle cuts, with the latter depending on the 
(rescaled) photon energy cutoff zq. The numerical value we chose in this illustration is given 
by Zq = 0.68 and corresponds to Eq = 1.6 GeV. As seen, the main effect of the bremsstrahlung 
corrections is to shift the function for i = NL (red solid line) and i = NV (black dotted curve) 
down by a factor of about 1.7 and 1.3, respectively. Also the aforementioned logarithmic 
singularity for i = NV can be observed for p ^ 0. There is no shift for i = NH (blue dashed 
line) since we set p equal to nib in our illustration (see 02.131) ). 

We remark that other values for the renormalization scale p than nif, lead merely to a 
shift of the three curves plotted in the left frame of figure [2] by the same amount up or down. 
The reason for this is that the three quantities given in flA.4p - flA.6l) have exactly the same 



/i- dependence. The same comment is also true for the three curves shown in the right frame 
of figure as can be seen from fl2.9l) - fl2.14l) . However, it is clear that a variation of p will 
change the relative importance with which each individual contribution enters the function 
G-j^{Eq, p). For example, the choice p = 1.2 m;, leads to Y'^'^'^^\zq, p) ~ 0, and hence the 
fermionic corrections will be dominated by the two functions with i = NL and i = NV. On 
the other hand, for smaller values of p than nib the situation can be reversed. 

Next, we compare the fermionic corrections at NNLO with the NLO result. Using the 
numerical values as{nib) = 0.22, = 3, = 1 and Ny = 1 (for the other input parameters 
we use the same values as before), we find 

G7s{Eo, nib) = 0.086 - 0.009 + . . . = 0.077 + . . . , (3.1) 

where the two numbers after the first equality sign correspond to the 0{as) and the fermionic 
0{al) contributions given in fl2.4p . and the ellipses denote the still unknown 0{al) terms 
proportional to Cp and Ca as well as higher order corrections. Thus, at p = nib, the effect of 
the NNLO fermionic corrections is to lower the NLO value of GjsiEo, nib) by around 10%. For 
p = 2.5 and 7.5 GeV, the 0{as) term in Gj^^Eq, p) changes to 0.009 and 0.124, respectively, 
and the 0{a1) correction shifts these values by around 3% and -11%, respectively. 

Finally, we estimate the effect of the fermionic corrections at 0{a1) on the branching 
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ratio oi B ^ Xg'j. As seen in figure HI the contributions with massive charm quark loops 
(i = NV) at the physical value p ~ (0.262)^ are of comparable size as those with massless 
quarks in the loops {i = NL). Hence, the charm quark mass effects can with good accuracy be 
described by a single massless quark entering the large- /3o approximation. On the other hand, 
approximating the contributions due to massive bottom quark loops {i = NH) by massless 
ones is not very accurate (see figure [2]). Here, however, one should bear in mind that in the 
physical application we have three massless and only two massive quarks. That is the leading 
correction to Br (5 X^'y) will be given by the sum of the contributions with i = and 
i = Nv, where the former is weighted by a factor of three, and the correction due to the 
contribution with i = Nfj will only appear at the subleading level. In fact, the exact result of 
the fermionic corrections at /i = can be accurately approximated by setting Nl = 3.6 and 
Nh = Nv = in (12.51) . Thus, the effect of the massive bottom quark loops at p = nib can be 
accounted for in the large-/5o approximation by reducing the number of massless quark flavors 
by -0.4. Given that the large-/3o corrections of the [Oj, 08)-interference affect the branching 
ratio of i? — > Xg'j by around 0.7% for p = mb, 'we conclude that this will not be altered 
drastically when implementing the exact results for the fermionic corrections with massive 
quarks0 We remark here that for other values of the renormalization scale fj, than nib it 
happens that the exact result can only be approximated by the massless contribution when 
using a negative number of massless quark flavors. For example, for /i = mb/2, the exact 
result can be approximated by setting A'^^: = —1-32 and Nh = Ny = in (12.51) . Determining 
the effect of the new fermionic corrections on Br(5 — >■ Xs'^) in this case would require to 
repeat the interpolation procedure performed in [6]. Since we expect that it will also be 1% 
at most (when combining large- /3o and massive quark loop corrections), that is below the 
total uncertainty given in (II. 2p . we postpone this to a forthcoming analysis which will also 
include other contributions not considered so far. 



4 Summary 

In this paper we calculated the NNLO fermionic corrections to the total decay rate and 
the photon energy spectrum induced by the interference of the electro- and chromomagnetic 
dipole operators. We confirmed the results for the 0(a^/?o) terms given in [14,16] and also 
presented analytical results for the contributions with massive bottom and charm quark 
loops. We expect that the combination of both the massless and the massive quark loop 
contributions affects Br(_B -^s7) by 1% at most. 
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A Two-particle cuts 

In this appendix we specify the contributions of the 2-particle cuts entering the functions 
r(^) and ^ith i = NL, NH, NV. To this end we write 

F«(zo,/x) = ril^lO,/.) + r3-cuts(^o,/x) (A.l) 
and 

y(2)(zo,/i) = 1^2?uts(0,/i) + llcL(^0,^) + 4cL(^0,/i) , (A.2) 

with F„^_cuts ^iid l^i.^cuts summing all contributions of the n-particle cuts. The contribution of 
the 2-particle cuts at 0{as) reads 

^icL(0, /i) = ^ (33 - 2rr') + ^ , (A.3) 
and those at 0{a'^) are given by 

.(2,nl) 



yr;;:::(0, Z^) = (157 - Sn^) - ^ (47 - 2n^) - f + f Ca , (A.4) 
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YlZ:\0, /^) = ^ (244 - 27v^vr - 61vr^) " ^ (47 - 2vr^) L, 

-y^J-gC3 + 32v/3Cl.(^), (A.5) 
y£;2\0, /^) = [157 - 279p - (5 + 9p^ - 42^^/^)] - ^ (47 - 2n^) L, 

- y + ^ (19 + 20p) Inp + ^ (8 + 9p') \n'p 

- ^ In^ - ^ (1 - 6/) $i(p) - ^(19 - 46p) Up) 

- p(13 + 14p) Up) - ^ Up) - f lnpLi2(l - p) 

+ I (5 + 9p^ + 14p^/^) Li,(l - p) - ip^p3/2 Li, (1 - ^p) 

64 64 / 1\ 

+ -Li3(l-p) + -Li3 1-- . (A.6) 



9 ^ " 9 V P 

In writing our results at O(a^) we tacitly performed a splitting of ^2-cuts(0; /^) i^^^ terms 
being proportional to certain combinations of colorfactors, in complete analogy to what we 
did in (12. 5p . We remark that all contributions given above are by themselves free of infrared 
and coUinear singularities, as well as independent of the gauge parameter entering the gluon 
propagator. 
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B Auxiliary functions 



Here we collect the four auxiliary functions $ri(p) introduced in section [2l They are defined 
as follows, 



$i(p) = 9{1 - 4p) [In^ y-n^]- e{Ap - 1) arccos^ ~ ^) ' 
$2(p) = Vll -4p| 1^(1 - 4p) Iny - ^(4p - 1) arccos - Yp 
$3(p) = V^T^(^(l-4p) 



Li2(-t/) + -ln'l/ + — 



1 



^(4p- l)Cl2 2arcsin 



2v^ 



1 3 ■ 
Li3 {-y) + Y^lii y + Y^lny 



$4(p) = ^(l-4p) 

+ ^(4p- 1)013 [2 arcsin 
where ^(2;) is Heavyside's step function 



2v^ 



y 



1 - yi - 4p 

1 + - 4p 



(B.l) 
(B.2) 



(B.3) 



:b.4) 



:b.5) 



and p > 0. The definitions of the two Clausen functions appearing in the above given 
equations read [32] 



CUz) = Im[Li2 (e*^)] , Ch{z) = RefLis (e"^^)] 
and those of the polylogarithms are given by 



- ln(l-a;) 
112(2;) = -/ dx , 113(2;) 



X 



X 



:b.6) 



(B.7) 
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